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Abstract 

We present a trivial probabilistic illustration for representation of 
quantum mechanics as an algorithm for approximative calculation of 
averages. 

1 Introduction 

In a series of authors papers there was developed an asymptotic ap- 

proach to the problem of hidden variables in quantum mechanics. It was 
shown that, in spite of all "NO-GO" theorems (e.g., von Neumann Cohen- 
Specker, Bell it is possible to construct a classical physical model, Pre- 
quantum Classical Statistical Field Theory - PCSFT, such that quantum 
mechanics can be considered as an approximation of PCSFT with respect to 
a small parameter a. This parameter is given by the dispersion of fluctuations 
of classical random field which is represented in quantum formalism by a von 
Neumann density operator. In PCSFT the role of hidden variables is played 
by classical fields. This is a model of classical statistical mechanics with 
the infinite dimensional phase space. The infinite dimension induces a num- 
ber of rather technical mathematical problems, in particular, using theory of 
measures on Hilbert spaces. Such purely mathematical difficulties embarrass 
understanding of PCSFT. In this paper there will be presented a simplified 
version of the approximation algorithm which induces the quantum rule for 
calculation of averages - the von Neumann trace-formula 0. 
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2 Method of Taylor approximation for func- 
tions of random variables 



Here we follow chapter 11 of the book [7| of Elena Ventzel. This book was 
the basic book for teaching probability theory in Soviet military colleges.^ 
Elena Ventzel wrote her book in the form of precise instructions what student 
should follow to solve a problem: 

"In practice we have very often situations in that, although investigated 
function of random arguments is not strictly linear, but it differs practically 
so negligibly from a linear function that it can be approximately considered as 
linear. This is a consequence of the fact that in many problems fluctuations of 
random variables play the role of small deviations from the basic law. Since 
such deviations are relatively small, functions which are not linear in the 
whole range of variation of their arguments are almost linear in a restricted 
range of their random changes," 0, p. 238. 

Let y = f{x). Here in general / is not linear, but it does not differ so 
much from linear on some interval [rrix — d,mx + S], where x = x{uj) is a 
random variable and 

is its average. Here (5 > is sufficiently small. Student of a military college 
should approximate / by using the first order Taylor expansion at the point 
: 

y{uj) ^ f{m^) + f\mx){x{uj) - m^). (1) 
By taking the average of both sides he obtains: 

my^fim^). (2) 

The crucial point is that the linear term f'{mx){x{u>) — rrix) does not give 
any contribution! Further Elena Ventzel pointed out p. 245: "For some 
problems the above linearization procedure may be unjustified, because the 
method of linearization may be not produce a sufficiently good approxima- 
tion. In such cases to test the applicability of the linearization method and 

am thankful to my farther-in-law, Alexander Choustov (marine artillery ofRcer) who 
pointed out to this chapter. 
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to improve results there can be applied the method which is based on pre- 
serving not only the linear term in the expansion of function, but also some 
terms of higher orders." 

Let y = f{x). Student now should preserve the first three terms in the 
expansion of / into the Taylor series at the point : 

y{uj) f{m^) + f'{m^){x{u) - m^) + ]^f" {mx){x{u) - m^f . (3) 

Hence 

my^f{m^) + ^f"{m^), (4) 

where 

a1 = E [x — m^Y = J (x^uj) — nij-Y (iP(a') 

is the dispersion of the random variable x. 

Let us now consider the special case of symmetric fluctuations: 

nix = 

and let us restrict considerations to functions / such that 

/(O) = 0. 

Then we obtain the following special form of Q: 

my ^ ^r(O). (5) 

We emphasize again that the first derivative does not give any contribution 
into the average. 

Thus at the some level of approximation we can calculate averages not by 
using the Lebesgue integral (as we do in classical probability theory), but by 
finding the second derivative. Such a "calculus of probability" would match 
well with experiment. I hope that reader has already found analogy with the 
quantum calculus of probabilities. But for a better expression of this analogy 
we shall consider the multi-dimensional case. Let now 

X , . . . , , 

SO we consider a system of n random variables. We consider the vector 
average: 

rrix = {rrix^, ...,m^„) 



3 



and the covariance matrix: 

Wc now consider the random variable y{uj) = f(xi{uj),...,Xn{uj)). By using 
the Taylor expansion we would like to obtain an algorithm for approximation 
of the average niy. We start directly from the second order Taylor expansion: 

1=1 

1 " d"^ f 

+ 2 5Z [m^,. ■■■,rn:,J){xi{u) - m:,^){xj{u) - m^.), (6) 

and hence: 



1 " 5^/ 

i,j=l * ^ 

By using the vector notations we can rewrite the previous formulas as: 

y{uj) f{m^) + {f{ma;),x{u)-m^) + ^{f'\m^){x{uj)-m^),x{u;)-m^). (8) 
and 

mj,«/(m^) + -TrB,f"(m,). (9) 
Let us again consider the special case: rrix — and /(O) = 0. We have: 

my ^ ^Tr Bxf"(0). (10) 

We now remark that the Hessian /"(O) is always a symmetric operator. Let 
us now represent / by its second derivative at zero: 

Then we see that, at some level of approximation, instead of operation with 
Lebesgue integrals, one can use linear algebra: 

my p^Tr BxA (11) 
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